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We discuss few selected topics related to the calculation of hadronic amplitudes 
o ■ relevant for two-body non-leptonic B decays. 

■ liitroductioii 

It is likely that most of the future results in B phenomenology, including the study 
of CP violation, will come from the measurements of two-body non-leptonic decays. 
The major theoretical problem in predicting the rates of these decays is the evalua- 
tion of hadronic amplitudes. In general the solution to this problem is unknown as 
it contains all the difficulties of low-energy strong interactions and hadronization. 
In QCD, the best one can do is using the operator product expansion to separate 
the short- and long-distance scales. In the resulting effective Hamiltonian, the effect 
of short-distance physics can be computed perturbatively and it is described by the 
notorious Wilson coefficients. Long-distance, non-perturbative physics is contained 



*Talk given by M.C. at Beauty '97, 5th Intl. Workshop on B Physics at Hadron Machines, 
UCLA, Los Angeles, CA, USA, October 13-17, 1997. 



1 



in the hadronic matrix elements of a set of local operators. In particular, for B 
decays into two mesons, one needs to compute matrix elements of dimension-six 
four-fermion operators (we neglect in the following magnetic dipole transitions), 

(M l M 2 \Q l \B) = (M l M 2 \bT l q l q 2 T' i q z \B) , i = 1, . . . , 10 , (1) 

where b, qj are the appropriate quark fields and Tj, 1^ are various combinations of 
Dirac and colour matrices. Unlike K physics, neither analytic computation tech- 
niques nor systematic expansions are known for these matrix elements. Numerical 
approaches to QCD are severely limited by the present computing power, since a 
small lattice spacing is required to simulate the heavy B field. 

These problems have not prevented theorists from making predictions of inter- 
esting quantities, such as BRs, asymmetries, etc. On the one hand, observables free 
from hadronic uncertainties have been identified, the best (and apparently unique) 
example being the CP asymmetry in Bj — > J/^Ks- On the other hand, when the 
evaluation of hadronic matrix elements could not be avoided, as for the BRs, vari- 
ous theoretical approaches, such as flavour symmetry, factorization and form-factor 
models, have been developed and used. 

Although these approaches have been successful in some applications, their the- 
oretical soundness is questionable. From a phenomenological point of view, this 
means that the theoretical error affecting the predicted amplitudes can only be 
guessed. 

Factorization 

To be concrete, let us consider the popular approach based on the factorization 
hypothesis applied to an emission-dominated decay, say B + — > D + n°. Using Fierz 
and colour rearrangement, one is left only with the matrix element 

(L> + 7T°|&7„(1 - cY(l - l 5 )d\B + ) , (2) 

where the quark fields are contracted according to the disconnected emission diagram 
DE in fig. [I]. The factorization hypothesis states that 

(D + n°\h,(l - l5 )u crf{l - -f 5 )d\B + ) = - ^)u\B + ) (D + \c^{l - 75 )d|0) , 

(3) 

namely the four-fermion operator matrix element is given by the product of the 
matrix elements of two currents. Pictorially this means that the two disconnected 
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Figure 1: Diagrams representing the Wick contraction of a four-fermion operator be- 
tween a B and two mesons Mi and M 2 . The operator insertion is represented by the 
double black dots. Disconnected emission (DE), penguin (DP) and annihilation (DA) 
diagrams are shown. 

branches of the emission diagram do not interact. The heuristic physical argument is 
appealing |1]: if a large energy is transferred to the emitted meson, as in the case at 
hand, the emitted quarks have not "enough" time to interact before going far from 
the interaction point and hadronizing. Starting from this picture, Dugan and Grin- 
stein introduced the large energy effective theory (LEET) H, in which factorization 
holds at the lowest order of a systematic expansion in powers of Kqcd/ E errii t. Besides 
the objections that can be raised even on the lowest-order result M, higher order cor- 
rections in LEET, or equivalently corrections to factorization, are not known. This 
is already a crucial point: even accepting that factorization holds in some effective 
expansion of QCD, there is no indication that it is an accurate approximation, and 
corrections at the level of 10-20%, or even more, are not surprising at all. 

Charming penguins 

It is easy to find decays where a moderate deviation of the matrix elements from 
their factorized values changes the predicted BRs by orders of magnitude. The 
reason is simply that some matrix elements, which vanish in the factorization limit, 
are Cabibbo-enhanced with respect to emission diagrams. This is precisely the 
mechanism which produces the large charming-penguin enhancement j|, ||. 

In general, it is useful to identify two classes of operators: i) current-current 
operators of the form bTq u q' u T's which have 0(1) Wilson coefficients and ii) pen- 
guin operators like bTsJ2 q <lF'q, the Wilson coefficients of which are ^ 0.03 at a 
scale ~ Mb- In B — > Kir and B —>■ Kp, for example, the only operators of class 
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Figure 2: Dependence of BR(B + — > p + K°) on t]l, the ratio of charming penguin to 
emission amplitudes. Exact factorization prediction corresponds to r] L = 0. 

i) which have non-vanishing factorized matrix elements are those containing up- 
quark fields. They are doubly Cabibbo-suppressed with respect to the operators in 

ii) , so that their Wilson coefficient enhancement is compensated and the contribu- 
tion of the two classes is comparable. However, if we allow for a violation of the 
factorization, class i) operators containing charm-quark fields can also have non- 
vanishing matrix elements. Contractions of these operators, like DP in fig. [I], are 
called charming penguins and contribute to the decays we are considering. Since 
they have large Wilson coefficients and are not Cabibbo-suppressed, their contribu- 
tion easily becomes dominant, already assuming corrections to factorization at the 
level of 10 — 20% Q. Striking examples of charming-penguin dominance are given 
by B — > pK channels. For instance, BR{B + — > p + K°), shown in fig. |2|, changes by 
three orders of magnitude for modest values of r/i, which is the ratio of charming- 
penguin to emission amplitudes. Exact factorization prediction, namely t]l = 0, is 
definitely unreliable for this decay. Notice that the value of the charming-penguin 
parameters suggested by B — > Kn measurements brings the BR(B — > Kp) near to 
their present experimental bounds ||. 

This discussion leads to the following conclusion: there is a class of rare B decays 
for which the use of factorized amplitudes gives very unstable predictions, which are 
drastically changed by moderate corrections to exact factorization. These decays, 
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including B — ► Kit/Kp/Kr], etc. have been extensively studied in ref. ||. 

The recent CLEO measurements of BR(B — > Kir) || actually call for some 
enhancement over the predictions obtained with factorized amplitudes. If charm- 
ing penguins have to explain the data, their value must be about 20-30% of the 
corresponding emission diagrams ||. Notice that the recent analysis of ref. 0, 
which claims to be able to reproduce the data using factorized amplitudes, some- 
how takes into account charm-loop effects in perturbation theory by enhancing the 
penguin-operator Wilson coefficients. The underlying physical process providing the 
enhancement is the same in the two approaches, but we believe that the perturbative 
treatment is not appropriate. 

A somewhat related argument, which has recently become popular, is the effect 
of the final state interaction (FSI) in B decays. In two-body decays, FSI is neglected 
in factorized amplitudes by definition. Also in this case, there are amplitudes that 
are neglected on the basis of factorization. For example, the annihilation diagrams, 
like DA in fig. [I], are usually neglected on the basis of the following argument: 
the annihilating quarks must be near enough for the weak current to annihilate 
them, implying a suppression proportional to the B wave function at the origin, 
namely a factor /b/Mb- However this and other diagrammatic arguments may 
not hold in presence of FSI, which mixes up different classes of diagrams Since 
recent theoretical estimates suggest, at variance with factorization, that the final- 
state interaction cannot be neglected in B decays, many phenomenological analyses 
relying on neglecting this or that amplitude on the basis of factorization should be 
reconsidered. FSI effects in B — > tttt, B — > K% have been studied in the analysis of 
ref. . Moreover, it has been recently shown [|H]] that rescattering effects invalidate 
the Fleischer-Mannel bound JTTJ] on cos(7) and affect bounds on new physics fll2"| . 

And all that 

Charming-penguin enhancement is not always effective: there are decay channels 
for which penguins are not present or not enhanced. Table [I] contains a list of 
measured emission-dominated channels. In these cases, factorized amplitudes are 
expected to give more reliable predictions, yet further assumptions are required. Us- 
ing Lorentz invariance, matrix elements of currents can be parameterized in terms 
of form factors, which however are known only in some special cases. If the external 
states are both heavy, HQET helps to express them in terms of few known quanti- 
ties, such as the heavy masses, the Isgur-Wise function slope p 2 , etc fl3fl . On the 
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Channel 


QCDSR 


NRSX 


Experiment 


BR x 10 5 


( p 2 = o.65) 


(P 2 = 1.1) 




B d -» 7T+D- 


301 


331 


310 ± 44 


B d -» vr+Z)*- 


323 


308 


280 ± 41 




794 


866 


840 ± 175 


5 d - p+D*- 


994 


949 


730 ± 153 




508 


534 


500 ± 54 


S+ -> ^+£>*° 


605 


567 


520 ± 82 




1015 


112 


1370 ± 187 


S+ -» p+L>*° 


1396 


1339 


1510 ±301 


S d -> iC° J/^ 


81 


75 


85 ± 14 


£ d K*°J/$> 


164 


173 


132 ± 24 


b+ _> j/* 


84 


78 


102 ±11 


S+ -» K* + J/y 


171 


180 


141 ± 33 


e 


0.47 


0.42 






0.42 


0.40 




X 2 /dof 


1.0 


1.0 





Table 1: Predictions for measured emissic 
NRSX and QCDSR, are considered for 
parameters are chosen according to the 
fitted parameters £ and 5^ are shown i 
xVdof. 



>n-dominated decays. Two form-factor models, 
different values of p 2 . Other required input 
central values of ref. The values of the 
n the two cases together with the values of 



other hand, heavy-light form factors need phenomenological models to be evaluated. 
Input-parameter dependence in the heavy-heavy case, and model dependence in the 
heavy-light one, introduce further theoretical errors, which sum up with the uncer- 
tainty on factorization and should be taken into account in the phenomenological 
analyses. 

Once form factors are chosen, no other free parameters are present in exact 
factorization and emission-dominated decay rates should be predicted without fur- 
ther assumptions. However it has become a standard procedure, starting with the 
well-known a\ and ci2 of BSW to introduce more parameters to account for 
possible deviation from factorization and fit them to the experimental data, see e.g. 
refs. P, |6], [15] . The comparison between these fitted parameters and the factorization 
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Figure 3: Dependence of x 2 /dof of the global fit to emission-dominated decays on the 
value of the Isgur-Wise function slope p 2 appearing in the heavy-heavy form factors for 
two different form-factor models. 



expectation provides a test of the reliability of factorization. Our parameterization 
is the following: connected (CE) and disconnected (DE) emissions are related ac- 
cording to CE = £e zSi DE, where £ and 5^ are the parameters to be fitted, then DE 
is expressed in terms of form factors using factorization. Colour rearrangement and 
exact factorization would give £ = 1/3 and 5% = 0. Fitting the decay channels in 
tab. p], we find, in agreement with other analyses f| [15], that experimental data are 
reasonably well described by factorized amplitudes for reasonable values of £ and 5%. 
Surprisingly enough, we have found a strong dependence on the Isgur-Wise function 
slope p 2 appearing in heavy- heavy form factors. Indeed, by changing p 2 , the value 
of x 2 m the fit also changes as shown in fig. |3|. In the global fit, the details of the 
model used for the heavy-light form factors turn out to be hidden by this large de- 
pendence. We have considered two popular models: NRSX fl6 |, which is a modern 
version of the original BSW model, and QCDSR |fl7j] , a model based on light-cone 
QCD sum rules. We have found that NRSX, QCDSR and the other models consid- 
ered in ref. || fit well the data, once a suitable value of p 2 is assumed. However, for 
example, the best fit within QCDSR and NRSX is given by quite different values of 
p 2 (p 2 ~ 0.65 and ~ 1.1 respectively). Moreover there exist experimentally allowed 
values of p 2 for which no model gives a good fit. 
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It is reassuring that it is always possible to fit the data with acceptable val- 
ues of the parameters. Within a given class of two-body final states with definite 
Lorentz properties, this may simply be a consequence of the dominance of emission 
diagrams and of SU(3) flavour symmetry, which reduce the number of independent 
amplitudes, rather than a test of the factorization hypothesis. However the global 
fit of many decay channels, including pseudoscalar and/or vector mesons in the final 
state, actually probes the different form-factor models and the factorization itself. 
On the other hand, the strong p 2 dependence in the fit calls for a careful study of the 
uncertainties affecting the factorized amplitudes. Indeed, the common procedure of 
comparing different heavy-light form-factor models for fixed value of p 2 , certainly 
leads to underestimate the theoretical error on £ and 5^ (or whatever parameters 
are fitted). On the contrary, p 2 should be varied within its experimentally allowed 
range in order to estimate the actual theoretical uncertainty. 
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